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ABSTRACT 


Two  techniques  were  developed  for  determining  the  kinetic  parameters  of  chemical 
reactions  from  differential  scanning  calorimeter  data.  In  one  technique,  variation  with  time 
of  the  heat  gain  or  loss  as  reactions  proceed  can  be  used  to  determine  the  order,  rate 
constant,  activation  energy,  and  heat  of  reaction.  In  the  second  technique,  variation  with 
time  of  the  glass  transition  temperature  of  certain  materials  can  be  used  to  determine  t  e 
order,  rate  constant,  and  activation  energy.  An  application  has  been  made  to  the  curing  of 
an  epoxy  resin. 
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SUMMARY 

Differential  scanning  calorimeters  (DSC)  measure  the  rate  at  which  energy  must  be  applied 
to  a  sample  or  to  a  reference  holder  to  maintain  them  at  the  same  temperature.  Two 
techniques  were  developed  for  obtaining  kinetic  data  using  a  DSC. 

If,  at  a  given  temperature,  a  sample  is  undergoing  a  reaction  involving  an  energy  change,  the 
output  of  the  DSC  is  proportional  to  the  rate  of  heat  flow  due  to  that  reaction.  Equations 
were  derived  in  which  the  variation  with  time  of  the  DSC  output  could  be  used  to  determine 
reaction  orders,  enthalpy,  rates,  and  activation  energies.  The  exothermic  curing  of  a  B-staged 
epoxy  resin  (ERL-2258/ZZL-0820)  was  studied  as  an  example  and  was  found  to  follow 
second-order  kinetics  with  a  rate  constant  of: 

k2  =  QBQe^^^'^^^^gcaf^  sec"^  . 

For  amorphous  materials,  a  change  in  heat  capacity  occurs  as  they  go  from  a  glassy  to  a 
liquid  state  at  the  glass  transition  temperature  (Tg).  The  Tg  of  some  materials  varies  in  a 
regular  manner  with  the  extent  of  polymerization.  Equations  were  derived  that  relate  this 
variation  of  Tg  with  the  kinetic  parameters.  Curing  of  the  B-staged  epoxy  resin  was  studied 
by  determining  the  Tg  from  DSC  heat-capacity  curves  for  samples  heat  treated  for  various 
times  and  temperatures.  The  results  were  found  to  correlate  vvell  with  those  obtained  using 
the  first  technique. 

The  extent  of  cure  was  calculated  using  kinetic  data  for  a  typical  slow  curing  cycle.  The 
extent  of  cure  calculated  from  the  Tg  measured  for  a  sample  heat  treated  over  that  cycle 
corresponded  closely  with  that  calculated  from  kinetic  data. 
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INTRODUCTION 

This  study  was  undertaken  at  the  Oak  Ridge  Y-12  Plant(a)  to  develop  techniques  which 
could  be  used  to  determine  the  kinetic  parameters  of  chemical  reactions  using  a  DSC. 

Exothermic  and  endothermic  reactions  involve  heat  transfers  at  rates  which  decrease  with 
time  in  the  same  proportion  that  the  concentrations  of  the  reacting  species  decrease.  The 
DSC,  with  its  fast  response  and  high  sensitivity,  is  well  suited  to  making  kinetic  studies  of 
such  reactions.  Previous  workersH)  have  reported  techniques  for  determining  kinetic 
parameters  from  data  obtained  under  nonisothermal  conditions.  Interpretation  of  the  results 
is  not  straightforward (2)  and  becomes  virtually  impossible  when  simultaneous  reactions 
with  different  activation  energies  occur.  Rogers(3)  found  that  plots  of  log  h  vs  t  gave 
straight  lines  for  first-order,  isothermal  reactions,  where  h  is  the  DSC  output  and  t  is  time. 
Gray(4)  determined  the  fraction  of  heat  evolved  at  t  by  planimeter  measurements  and 
equated  it  with  the  fraction  of  cure  of  an  epoxy  resin.  Comparison  was  then  made  with  the 
Tg.  In  this  work,  equations  were  derived  which  relate  h,  taken  under  isothermal  conditions, 
with  the  kinetic  parameters.  They  were  applied  in  the  determination  of  rate,  order, 
activation  energy,  and  heat  of  the  curing  reaction  of  an  epoxy  resin. 

For  many  polymers,  Tg  varies  in  a  regular  manner  with  molecpiar  weight. (5)  Equations  were 
derived  which  relate  Tg  changes  with  the  kinetic  parameters  and  were  used  to  determine  the 
rate,  order,  and  activation  energy  of  the  curing  of  an  epoxy  resin. 

The  utility  of  the  data  reported  herein  is  that  they  may  be  used  to  calculate  rather  than 
measure  the  degree  of  cure  obtained  from  a  variety  of  proposed  heat-treatment  cycles.  An 
example  of  such  a  curing  cycle  is  given  in  which  the  calculated  degree  of  cure  was  found  to 
correspond  closely  with  the  measured  degree. 

Effects  of  changes  of  stoichiometry  and  addition  of  inert  diluents  upon  the  determination 
of  the  kinetic  parameters  were  also  studied. 


(a)  Operated  by  the  Union  Carbide  Corporation's  Nuclear  Division  for  the  US  Energy 
Research  and  Development  Administration. 
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KINETICS  OF  CURING  AND  HEAT  CAPACITY  OF  AN  EPOXY  RESIN 
APPARATUS  AND  TECHNIQUES 

Thermal  measurements  and  heat  treatments  were  performed  on  a  Perkin-Elmer  Model 
DSC-2  differential  scanning  calorimeter  (DSC).  This  instrument  contains  a  matched  pair  of 
9-mm-diameter  sample  holders  with  integral  heating  elements  and  platinum  resistance 
thermometers.  Heating-element  power  is  derived  alternately  at  85  Hz  from  two  sources.  The 
average  temperature  of  the  two  holders  is  controlled  by  a  feedback  loop  and  a  digital 
temperature  programmer  controlling  one  power  supply.  The  temperature  difference 
between  the  holders  is  nulled  by  the  application  of  power  from  the  second  supply  and 
feedback  loop  to  the  holder  at  the  lower  temperature.  The  power  from  this  second  source, 
which  is  the  time  derivative  of  the  transition  or  reaction  energy  or  of  the  heat  capacity  of 
the  sample,  is  recorded  as  a  function  of  time  or  temperature.  Samples  of  0.1  to  200  mg  can 
be  studied  over  the  range  from  100  to  1000  K  at  heating  or  cooling  rates  from  0  to 
320  K  min”1  and  at  sensitivities  down  to  0.01  meal  sec"1 .  Instrument  response  is  faster,  t  e 
sample  size  smaller,  and  the  output  more  physically  meaningful  than  with  conventional 
differential  thermal  analysis  instruments. 

In  this  work,  sample  sizes  ranged  from  4  to  20  mg.  Glass  transition  temperatures  (Tg)  were 
determined  at  heating  rates  of  5  deg  min'l  and  were  defined  as  the  temperature  at  which 
one  half  the  change  in  heat  capacity  (difference  between  the  glass  and  liquid  states)  had 
occurred.  The  heat  of  fusion  of  indium (6)  at  429.8  K  and  the  heat  capacity  of  aluminurn 
oxide(7)  were  used  as  calibrating  standards  for  the  determination  of  the  reaction  and 
transition  enargias  and  of  the  heat  capacity  of  the  resin. 

The  epoxy  resin  used  in  this  study  was  a  31:9  mixture  of  ERL-2258/ZZL-0820_(b) 
ERL-2258  consists  of  50%  diglycidyl  ether  of  bisphenol-A  (DGEBA)  and  W/o 
bis(2,3-epoxycyclopentyl)  ether.  ZZL-0820  consists  of  53%  m-phenylenedtamine,  34/o 
p,p'-methylenedianiline  and  13%  DGEBA.  The  mixture  was  heat  treated  at  358  K  for  two 
hours  to  effect  a  partial  cure,  or  B-stage. 

RESULTS  AND  DISCUSSION 

Three  kinds  of  data  were  obtained  from  the  DSC  thermograms:  (1 )  rates  of  heat  evolution 
from  the  curing  reaction  of  the  resin,  (2)  relative  and  absolute  heat  capacities,  and 
(3)  transition  energies.  Rates,  the  order,  and  the  activation  energy  of  the  curing  reaction 
were  calculated  from  the  rates  of  heat  evolution  at  different  temperatures.  Glass  transition 
temperatures  were  obtained  from  the  relative  heat  capacities. 

Determination  of  Kinetic  Parameters  by  Heat- Evolution  Methods 

Samples  were  heated  to  423  or  473  K  at  320  deg  min-l  -  and  the  ratethat  heat  evolved  per 
gram,  h,  from  the  isothermal  exothermic  curing  reaction  was  plotted  as  a  function  of  time. 


(b)  A  product  of  the  Union  Carbide  Corporation 
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Figure  1.  RATE  OF  HEAT  EVOLUTION  AT  TWO  TEMPERATURES  FOR  CURING  AN  EPOXY  RESIN.  (The  Initial 
Irregularity  in  the  Curve  of  the  473  K  Run  is  Due  to  Slight  Temperature  Changes  in  the  Heat  Sink) 


Figure  1  is  a  graph  of  the  typical  data  obtained.  During  the  initial  400  seconds  of  the  473  K 
run,  the  instrument  heat  sink  was  not  close  enough  to  thermal  equilibrium  for  accurate  data 
to  be  taken  at  the  high  sensitivities  used.  At  5000  seconds  into  the  423  K  run  and  at  3800 
seconds  into  the  473  K  run,  the  uncertainty  in  the  rate  of  heat  evolution  increased  because 
of  very  slight  instrument  drift,  and  the  runs  were  therefore  terminated  before  completion  of 
cure.  The  area  under  each  of  the  two  curves  should  be  equal  to  the  heat  of  the  curing 
reaction,  AH.  Planimeter  determinations  of  such  areas  were  found  to  vary  greatly,  largely 
because  of  the  inaccurate  curve  traced  at  the  start  of  each  run. 

However,  the  shape  of  the  major  portions  of  the  curves  can  be  analyzed  by  equations 
derived  in  Appendix  A  for  first,  second,  and  third-order  kinetics  to  yield  AH,  the  activation 
energy  E+,  and  the  rate  constants  for  the  curing  reaction,  kp.  For  first-order  reactions  plots  of 
In  h  vs  t  give  straight  lines  according  to  Equation  A-1 ;  but,  with  this  epoxy  resin,  such  plots 
gave  curved  lines,  indicating  that  the  curing  reaction  is  not  first  order.  Plots  of  h~1/2  vs  t 
give  straight  lines  for  second-order  reactions  according  to  Equation  A-7.  As  noted  in  Figure 
2,  these  data  lie  on  straight  lines,  demonstrating  that  second-order  kinetics  are  followed.  The 
slight  curvature  after  3500  seconds  for  the  423  K  heat  treatment  is  due  to  a  very  slight 
instrument  drift.  These  points  were  not  used  in  calculations  of  slope  and  intercept.  The 
points  taken  before  400  seconds  were  likewise  not  used  because  thermal  equilibrium  of  the 
heat  sink  with  the  sample  holders  had  not  yet  occurred. 
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Figure  2.  SECOND-ORDER  KINETIC  PLOT  OF  THE  RATE-OF-HEAT-EVOLUTION  DATA  AT  TWO  TEMPERA¬ 
TURES  OF  HEAT  TREATMENT.  (Data  Taken  at  Less  Than  500  Seconds  and  More  Than  4000  Seconds  were  not  Used  in 
Computing  the  Least-Squares  Lines) 


The  second-order  reaction  constant,  k2,  was  computed  from  the  squares  of  the  slopes  of  the 
lines  in  Figure  2  and  was  found  to  have  values  of  2.07  x  lO'S  and  1.34  x  10“ g  cal  1  sec  1 
at  423  and  473  K,  respectively.  The  activation  energy  is  thus  14.8  kcal  morL  and  the 
second-order  rate  constant  at  temperature  T  is  given  by: 


k2  =  959e~^'^'^^^^gcal  ^  sec  ^ 


(1) 


The  reaction  half  times  were  computed  from  the  ratio  of  the  intercept  to  the  slope  and  were 
1664  and  202.4  seconds  at  423  and  473  K,  respectively.  The  initial  rate  of  heat  evolution  is 
the  reciprocal  of  the  square  of  the  intercept  and  is  1 .742  x  10“2  qal  g“1  sec“l  at  423  K  and 
1.821  X  10“  1  cal  g“1  sec“l  at  473  K.  The  area  under  the  curve  of  Figure  1  is  AH,  which  is 
the  reciprocal  of  the  product  of  the  slope  and  intercept  of  the  lines  in  Figure  2,  and  was 
computed  to  be  29.0  cal  g“'' .  Planimeter  measurement  of  that  area  gave  21 .3  cal  g“1,  but 
the  kinetic-data  computation  is  more  reliable  because  of  the  uncertainties  in  h  during  the 
initial  part  of  the  runs. 


Using  the  method  of  Uricheck.CI )  attempts  were  made  to  obtain  kinetic  parameters  by 
measuring  the  areas  between  the  initial  and  fully  cured  heat-capacity  curves  obtained  by 
scanning  at  1.25  deg  min“''  over  the  range  from  323  to  523  K.  A  sensitivity  of  0.02  meal 
sec“  1  for  a  full-scale  recorder  deflection  was  required.  The  area  is  approximately  the  sum  of 
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the  heat  of  reaction  and  the  product  ACpAT,  where  ACp  is  the  difference  in  heat  capacity 
of  the  uncured  and  cured  resin  and  AT  is  the  difference  between  Tg^o  and  Tqq,  the  final  and 
initial  Tg  values.  However,  as  a  significant  portion  of  the  curing  reaction  had  not  taken  place 
during  this  heat  treatment,  as  the  slower  scan  speeds  required  to  get  complete  reaction 
require  even  greater  sensitivity,  and  as  higher  final  temperatures  could  not  be  used  without 
degrading  the  resin,  this  effort  to  obtain  AH  had  to  be  abandoned. 

Determination  of  Kinetic  Parameters  by  the  Glass-Transition-Temperature  Method 

Amorphous  materials  exhibit  discontinuous  changes  in  their  heat  capacities  and  thermal 
expansions  at  the  temperature  at  which  they  change  from  the  glassy  to  the  liquid  state.  This 
temperature  (actually  a  temperature  range)  is  called  the  "glass  transition  temperature",  Tg, 
here  defined  as  that  temperature  at  which  one  half  the  change  in  heat  capacity  occurs  during 
a  scan  of  increasing  temperature.  Because  Tg  is  slightly  dependent  on  the  rate  of 
temperature  increase  (due  to  finite  thermal  conductivity),  a  slow  scan  speed  of  five  degrees 
per  minute  was  used. 

The  heat  capacity  of  the  as-received  B-staged  resin  is  shown  in  Figure  3.  Heat  capacities  of 
the  B-staged  resin  above  373  K  cannot  be  directly  determined  because  the  heat  given  off  in 
the  curing  reaction  significantly  changes  the  apparent  baseline.  Figure  4  shows  the  heat 
capacity  after  a  heat  treatment  of  25  hours  at  423  K.  During  the  initial  scan  of  each  sample 
of  B-staged  resin  through  the  Tg,  a  large  endothermic  spike  was  always  found,  with  an  area 
of  0.838  cal  g"l.  Subsequent  scans  were  devoid  of  this  feature  and  exhibited  the  normal 
S-shaped  curve  from  which  the  Tg  could  be  determined.  This  phenomenon  is  probably  a 
result  of  a  bulk  stress  relaxation  of  the  material  or  a  molecular-level  free-volume  effect. (8) 
The  endothermic  spike  was  also  seen  only  on  the  initial  scan  through  the  glass  transition 
region  of  the  heat-treated  sample  in  Figure  4  Note  that  Tg  is  greater  than  the 
heat-treatment  temperature,  indicating  that  the  curing  reaction  continues  at  temperatures 
below  the  Tg,  here  accounting  for  an  increase  in  percent  reaction  (from  86  to  94%),  as 
determined  from  Tg  and  kinetic  data.  An  endothermic  spike  was  also  always  found 
whenever  the  material  was  held  slightly  below  the  Tg  for  a  significant  period.  Its  area 
appeared  to  increase  with  time.  Kinetic  data  for  the  curing  reaction  below  the  Tg  could  have 
been  determined  (vide  infra),  but  were  not,  due  to  the  lack  of  time.  It  is  estimated  that  the 
reaction  rates  are  decreased  one  or  two  orders  of  magnitude  when  the  heat-treatment 
temperature  is  less  than  the  Tg.  The  effects  of  this  continuing  cure  upon  other  physical 
properties  should  be  investigated  to  determine  the  suitability  of  use  of  this  material  in  each 
potential  application,  particularly  if  optimum  physical  properties  are  obtained  from  material 
cured  less  than  1 00%. 

The  Tgs  for  many  polymers  have  been  shown  to  be  a  linear  function  of  the  inverse  of  their 
molecular  weights.(5)  Using  this  assumption,  expressions  are  derived  in  Appendix  B  relating 
the  Tg  changes  with  nth-order  reaction-rate  constants.  For  the  resin  studied  here,  a  Tg  of 
473  K  was  the  highest  obtained  after  an  overnight  heat  treatment  at  473  K  (greater  than 
98%  cured)  and  was  thus  used  as  Tg^o-  Use  of  higher  heat-treatment  temperatures  for  long 
time  periods  showed  that  Tg  went  through  a  maximum  at  about  473  K,  indicating  that 
degradation  was  occurring.  Figure  5  is  a  first-order  plot  of  the  data  listed  in  Table  1  which 
was  obtained  by  finding  the  Tg  of  samples  heat  treated  at  373,  423,  or  473  K  for  various 
times.  The  high  curvature  of  the  lines  demonstrates  that  the  curing  reaction  is  not  first 


m  the  Second  Heating} 


of  the  Material  through 
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order.  Figure  6  is  a  second-order  plot  of  the 
data.  Here,  the  curvatures  are  small, 
indicating  second-order  kinetics  are  fol¬ 
lowed  and  are  due  to  the  uncertainties  in 
the  final  Tgoo.  The  slopesare  equal  to  k2d/m, 
where  k2  is  the  second-order  rate  constant, 
d  is  the  density,  and  m  is  the  slope  of  the  Tg 
vs  M"1  plot  (M,  mol  wt).  An  activation 
energy  of  1 1 .9  kcal/mol  was  calculated  from 
the  slope  of  the  line  in  Figure  7. 

Using  Equation  A-10  of  Appendix  A,  the 
reaction  half  time  at  373  K  was  calculated 
to  be  17,720  seconds.  Equation  A-8  was 
used  to  calculate  from  the  heat-evolution 
data  the  extent  of  reaction  for  the  samples 
listed  in  Table  1.  The  Tgs  were  plotted 
against  these  percentages  of  reaction  in 
Figure  8.  No  assumptions  concerning  the 
relationship  of  Tg  with  molecular  weight 
were  used  in  this  figure.  The  closeness  of  all 


Table  t 

VARIATION  OF  THE  GLASS  TRANSITION 

temperature  with  heat 

TREATMENT 


Heat  Treatment 


At  373  K  At  423  K  At  473  K 


Time 

(sec) 

Tg(1) 

(K) 

Time 

(sec) 

Tg 

(K) 

Time 

(sec) 

Tg 

(K) 

0* 

323.7 

0 

325.9 

0 

328.1 

75* 

324.4 

120 

336.8 

60 

355.7 

255* 

326.6 

420 

353.7 

120 

370.1 

735* 

332.6 

1,020 

372.7 

202 

377.6 

1,395* 

336.4 

1,620 

385.3 

300 

392.2 

1,969 

342.6 

2,495 

396.2 

472 

399.8 

3,180* 

349.5 

3,881 

406.3 

810 

413.9 

3,360 

348.6 

6,780 

416.3 

1,822 

431.5 

4,560 

357,7 

55,800 

441.7 

3,900 

447.1 

5,907 

361.8 

90,000 

444.2 

64,800 

473.0 

6,795* 

364.3 

61,200* 

393.2 

69,900 

393.2 

(1 )  Those  marked  with  an  asterisk  are  the  data  that  were  ob¬ 
tained  at  20  deg  min""! .  The  other  data  were  obtained  at 
5  deg  min~^. 


but  three  points  to  the  line  is  a  measure  of 
the  quality  of  the  kinetic  data.  Those  three  points,  marked  "a",  correspond  to  samples 
which  were  heat  treated  at  times  sufficiently  long  for  the  Tg  to  eventually  exceed  the 
heat-treatment  temperature  as  the  curing  reaction  proceeded.  The  unknown,  but  greatly 
decreased,  reaction  rate  resulting  from  this  phase  change  was  not  determined  in  this  work 
due  to  a  lack  of  time  and  was  thus  not  used  in  the  extent  of  the  reaction  calculations. 


Comparison  of  the  Two  Methods 

Equation  B-3  (Appendix  B)  was  used  to  calculate  the  extent  of  the  reaction  based  on  Tg 
data.  Figure  9  is  a  plot  of  the  extent  of  the  reaction  obtained  from  heat-evolution  data 
versus  that  obtained  from  Tg  data.  The  assumption  that  Tg  is  a  linear  function  of  the 
reciprocal  of  the  molecular  weight  is  not  valid  for  this  system  over  the  whole  interval  of  the 
curing  reaction,  as  the  points  in  Figure  9  do  not  lie  on  the  line  of  unit  slope.  The  straight 
portion  of  the  line  from  0  to  80%  reaction  indicates  that,  over  this  range,  Tg  is  proportional 
to  M"!.  If  it  is  assumed  that  Tgg  was  measured  incorrectly  and  that  a  Tgoo  value  of  ^^3  K 
does  not  apply  to  this  resin  system  at  less  than  80%  cure,  new  values  (TgQ  and  Tgoi,)  can  be 
calculated  from: 

Tqo  =  t  ITg....-Tgol,or  12) 


Tqq 


^y(ioo) 

+  :,oo" ' 'g«>“  'go' ' 


(3) 


where  y(0)  and  y(100)  are  the  end  points  of  the  straight  line  in  Figure  9  which  were 
determined  by  least-squares  analysis.  Using  the  computed  values  of  333.1  and  434.4  K  for 
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Figure  6  SECOND-ORDER  KINETIC  PLOT  OF  GLASS-TRANSITION-TEMPERATURE  DATA. 
(Slopes  of  the  Lines  are  k2d/m;  the  Slight  Curvature  of  the  Lines  is  Due  to  an  Inaccuracy  in  the 
Determination  of  Tg^) 


Tn '  and  Tn^.  respectively,  the  data  were  replotted  in  Figure  10  in  which  the  points  of  le^ 
tha°n  80%  reaction  now  lie  on  the  line  of  unit  slope.  Physical  interpretation  of  these  data 
without  other  data  is  not  possible,  but  it  does  suggest  that  there  are  5°^  Para l  ei 
competing  chemical  reactions  occurring  in  the  last  20%  of  the  cure  which  alter  the  structur 
of  the  material.  This  possibility  is  not  surprising  in  view  of  the  complex  nature  of  the 
components  of  the  resin. 

Figures  1 1  and  1 2  are  the  same  as  Figures  6  and  7  except  that  Tgo  and  Tg^  were  u^d  in  the 
calculations.  Note  that  the  lines  in  Figure  1 1  are  straighter  than  those  in  Figure  6.  The  s  ope 
of  Figure  12  was  used  to  calculate  an  activation  energy  of  14.4  kcal  mol  1,  m  excellent 
agreement  with  that  calculated  from  the  heat-evolution  kinetic  data. 

Application  of  the  Kinetic  Parameters  to  the  Heat-Treatment  Cycles 

In  Appendix  C  an  example  is  given  of  the  use  of  kinetic  data  to  detetmine  the  extent  of 
reaction  for  a  curing  heat  treatment  consisting  of  a  1.25  deg  min'l  incr^se  from  323  to 
423  K.  a  30-minute  hoid  at  423  K,  and.a  1.25  deg  mm  '  decrease  from  423  to  323  K.  A 
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Simpson's-rule  integration  of  Equation  C-2 
was  performed  for  the  heat  treatment  of  slow 
heating  to  423  K  and  holding  at  423  K  for  30 
minutes.  From  this  integration,  the  curing 
reaction  was  calculated  to  be  64.0%  com¬ 
plete;  which,  from  Figure  8,  corresponds  to  a 
Tg  of  399  K.  From  this  value  it  was  estimated 
that  the  Tg  would  be  about  403  K  after  the 
cooling  portion  of  the  heat  treatment.  A 
second  Simpson's-rule  integration  was  per¬ 
formed  using  403  K  as  the  lower  limit  on  the 
cooling  portion,  and  this  integration  resulted 
in  a  calculation  of  68.5%  completion  of  cure. 
The  calculations  were  checked  by  finding  the 
percent  cure  from  the  Tg.  During  the  deter¬ 
mination  of  the  Tg  of  the  sample  after  the 
heat-treatment  cycle,  a  large,  very  sharp, 
endothermic  spike  was  noted  in  the  Tg  re¬ 
gion.  A  second  scan  was  performed,  indica¬ 
ting  a  Tg  of  407  K,  which  corresponds  to  a 
73.6%  cure.  The  agreement  between  calcula¬ 
tions  is  very  good. 


higure  7.  DETERMINATION  OF  THE  ACTIVATION 
ENERGY  FROM  GLASS-TRANSITION-TEMPERATURE 
DATA.  (The  Activation  Energy  is  the  Negative  of  the 
Slope  Divided  by  the  Gas  Constant,  R,  and  is  11 .9  kcal 
mor’') 


Figure  8.  CHANGE  IN  THE  GLASS  TRANSITION  TEMPERATURE  WITH  THE  DEGREE  OF  CURE.  (The  Degree  of 
Cure  was  Calculated  from  Heat-Evolution  Kinetic  Data  and  the  Time  the  Samples  were  Heat  Treated  at  373,  423,  and 
473  K;  the  Points  Marked  "a"  Correspond  to  Samples  whose  TgS  are  Smaller  than  their  Heat-Treatment  Temperatures) 


L 


Percent  Cure  from  Modified  Gloss  Transition  Temperature  Data  3-  I  I!  Percent  Cure  from  Gloss  Transition  Temperature  Data 
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Time  (ksec) 


Time  (ksec) 

Figure  11.  SECOND-ORDER  KINETIC  PLOT  OF  REVISED  GLASS-TRANSITION- 
TEMPERATURE  DATA.  (The  Data  of  Figures  were  Replotted  with  Tg^  =  333.1  K  and 
Tg  =  434.1  K;  the  Slopes  are  k2d/m,  the  Curvatures  of  the  Data  Points  of  Figures  are  not 
Present  Here) 


Effects  of  Changes  of  Stoichiometry  and  Inert  Diluents 

During  the  preceding  discussion  it  was  assumed  that  the  B-staged  epoxy  resin  consisted  only 
of  an  epoxide,  E,  a  stoichiometric  amount  of  an  active  hydrogen  compound,  H,  and  some 
condensed  polymer  product,  P.  Because  the  results  were  obtained  on  a  per-gram  basis,  the 
effect  of  the  presence  of  P,  or  any  other  inert  additive  for  that  matter,  would  be  that  of 
a  diluent  (ie,  to  decrease  AH  by  a  factor  of  1/(1  +  F),  where  F  is  the  weight  fraction  of 
additive,  and  to  decrease  reaction  rates  by  a  factor  of  1/(1  +  F)!^,  where  n  is  the  order  of 
the  reaction).  Neither  the  determination  of  the  order  of  the  reaction  nor  of  the 
activation  energy  would  be  affected. 
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The  rate  of  polymer  formation  for  the 
second-order  kinetic  results  obtained  here  is 
d[P]/dtand  is  given  either  by: 

^  =  k2[E]  [H]  ,  (4) 

dt  ^ 

—  =k2[E]2,or  (5) 

dt 

—  =  k2[H]2,  (6) 

dt 


where  the  brackets  designate  the  concentra¬ 
tions  of  the  species.  The  rate  is  measured 
directly  by  the  DSC,  but  it  is  impossible  to 
distinguish  between  Equations  4,  5,  and  6  on 
the  basis  of  one  formulation  alone.  If  a 
weight,  AW^,  of  E  were  added  in  excess  of 
stoichiometric  requirements,  the  concentra¬ 
tion  of  E  would  then  be: 


[E] 


AWe 
1  + - 

_ Me _ 

Me  +  AWe  Mh  ' 

- +  — 

dE  dH 


(7) 


and  the  concentration  of  H  would  then  be: 


Figure  12.  DETERMINATION  OF  THE  ACTIVATION 
ENERGY  FROM  REVISED  GLASS-TRANSITION- 
TEMPERATURE  DATA.  (The  Activation  Energy  is  the 
Negative  of  the  Slope  Divided  by  the  Gas  Constant,  R, 
and  is  14.4  kcal  morl,  in  Good  Agreement  with 
Computations  Obtained  from  Heat-Evolution  Data) 


[H] 


1 

Me  +  AWe  Mh 
- — —  - 

dE  dH 


(8) 


where  Me  3nd  Mh  ^be  molecular  weights  and  dE  srid  dH  3*^6  the  densities  of  E  and  H, 
respectively. 

Rate  Expressions  4,  5,  and  6  tlien  become: 


d[P] 


=  k2 


1  + 


AWe 

Me 


[Me  +  AWe 

y  dE  “ 


dt 


(9) 
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d[P] 

dt 


k2 


[Mg  +  AW^ 

r^E 


2 


,  and 


dt 


=  k2 


1 

/Me  +  AWe  ^  Mh\^ 
\  dE  ^  dH/ 


(10) 


(11) 


respectively,  and  determination  of  the  rates  of  reaction  of  two  different  formulation  ratios 
of  E  and  H  would  demonstrate  which  rate  law  was  followed.  The  heat  of  reaction  per  gram, 
AH,  is  decreased  by  a  factor  expressed  as: 

1 


Me  +  awe  M|-i 
dE  dH 


The  effects  of  inert  diluents  on  the  Tg  of  polymers  cover  a  wide  range  from  no  effect  to 
large  changes,  which  generally  decrease  the  value  of  Tg.  Changing  the  ratio  of  [E]  to  [H] 
from  the  stoichiometric  ratio  would  decrease  the  difference  between  Tgg  and  Tgoo- 


I 
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APPENDIX  A 

DERIVATION  OF  EQUATIONS  FOR  DETERMINING  KINETIC  PARAMETERS  FROM 
ISOTHERMAL  DIFFERENTIAL  SCANNING  CALORIMETRY  DATA 

Figure  1  shows  examples  of  DSC  thermograms  for  exothermic  reactions;  here,  the  curing  of 
a  B-staged  epoxy  resin.  If  each  molecule  reacting  emits  or  absorbs  the  same  quantity  of 
energy,  then  the  area  under  the  curve  up  to  time  t  is  the  amount  of  reaction  which  has 
occurred  and  is  given  by: 

X  =  /*  f(h)dt , 

Jo 

where  h  is  the  instantaneous  heat  flow  rate.  The  total  area  under  the  curve  is  the  heat  of 
reaction,  AH; 


First-Order  Reactions 

For  first-order  reactions; 


AH 


f(h)dt  . 


k-|t  =  In 


f  — 1 

[ah- X J  ' 


where  ki  is  the  first-order  rate  constant.  If  it  is  assumed  that  f  (h)  is  of  the  form; 

h  =  hoe” then 


(A-1) 


X  = 


f  hoe”*^!^  dt  =—  (1  -  e”'^1^),and 

Jq  ki 


f 


AH  =  /  hne-^lt  dt  =  — 


ki 


Thus: 


k-j  t  =  In 


"o 

ki 


-  i!2  (1  -  e-kl^). 

.ki  ki 


-  In  -  =  ki  t 

e-klt  ^ 
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Therefore,  Equation  A-1  is  the  proper  form  for  first-order  reactions.  Plots  of  In  h  vs  t  give 
straight  lines  of  slope  -ki  and  intercept  In  ho  =  ln(AH  ki).  Rogers(3)  determined  ki  from 
In  h  vs  t  plots,  but  missed  the  opportunity  to  calculate  AH. 


Second-Order  Reactions 


For  second  order  reactions,  with  either  initial  concentrations  of  the  two  reacting  species 
equal  or  with  second-order  dependence  on  one  specie: 


k2t  = 


X 

AH  (AH  -  X)  ' 


{A-2) 


where  k2  is  the  second-order  constant.  If  f(h)  is  of  the  form: 

h  =  ho  (1  +gtr^  , 

where  g  is  a  parameter  to  be  determined,  then: 


X 


•'o 

/ 


_9  M 

(1  -Hgt)  ^dt  =  -—  and 
1  +gt 


AH  =  y  ho  (1  +gtr^dt  =  ^ 


(A-3) 


(A-4) 


Substituting  Equations  A-3  and  A-4  into  Equation  A-2  gives:' 


Thus: 


gt  g  (1  +gt)  ^  ^ 
1  -(-  gt  ho  ho 


1 

(1  +gt)2 


1 


•<2 


1 

-  +  t 

Lg 


,  and 


(A-5) 


(A-6) 


21 


Equation  A-5  can  be  rewritten: 


1 

h1/2 


■1/2  ^2 
=  k2  '  \  + 


1/2 


g 


(A-7) 


and  a  plot  of  h  vs  t  will  give  a  straight  line  of  slope  and  intercept 
The  fraction  reacted  after  time  t  is: 


X  gt  AH  l<2t 
AH  ~  1  +gt  ~  1  +  AH  k2t ' 


(A-8) 


from  which  it  is  easily  seen  that  g  is  the  reciprocal  of  the  half  time. 

From  the  Arrhenius  equation: 

kn(T)  =  ,  (A-9) 

where  E"*"  is  the  activation  energy  and  R  is  the  gas  constant,  the  temperature  dependence  of 
the  g  parameter  is  easily  shown  to  be: 

g(T)  =  goe  .  (A-10) 

- 1 

The  activation  energy  is  normally  found  from  the  slope  of  a  plot  of  In  k^  vs  T 


Third-Order  Reactions 

Assume  a  solution  of  the  form: 

h  =  ho  (1  +  Itf^^^  . 
where  I  is  a  parameter  to  be  determined. 


(1  +ltr^^^dt 


c 


AH 


(1  -utr^/^dt  = 


2K 


Thus: 
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For  third-order  reactions: 


k3t  =  1/2 


(AH  -  X)2  AH 


Then: 


k3  - - 5  '  or 

8ho2 


LskaJ 


(1  +  It)' 


Upon  raising  this  expression  to  the  -2/3  power: 

1/3 

^-2/3  -  (i+it)  =  2k3^/^t-r2k3^/^r‘'  . 

Plots  of  h"^^^  vs  t  give  straight  lines  of  slope  2k3^/2  and  intercept  2k3^/^l  .  The 

slope-intercept  ratio  is  I.  The  fraction  reacted  at  time  t  is. 

—  =  1  -  (1  +  it)"^^^ , 

AH 

from  which  it  can  be  shown  that  the  concentration  of  the  reacting  species,  c,  is: 


where  t^  /2  time  of  the  reaction. 

Similar  methods  can  be  used  to  derive  equations  for  other  reaction  orders.  It  should  be 
noted  that  uncertainties  in  the  time  of  initiation  will  not  affect  the  determination  of  the 
rate  and  activation  energy.  Corrections  of  the  other  parameters  can  be  made  easily  by 
computing  intercepts  from  values  of  AH  determined  by  standard  techniques. 
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APPENDIX  B 

DERIVATION  OF  EQUATIONS  FOR  DETERMINING  KINETIC  PARAMETERS  FROM 
GLASS-TRANSITION-TEMPERATURE  DATA 

The  glass  transition  temperature,  Tg,  is  that  temperature  (or,  more  correctly,  temperature 
range)  at  which,  upon  heating,  an  amorphous  material  changes  from  a  glassy  to  a  liquid 
state.  It  can  be  determined  from  heat-capacity  or  thermal-expansion  curves. 

For  many  polymers,  the  empirical  formula: 

Tg  =  mM~^  -t  b 

has  been  found,  (3)  where  M  is  the  molecular  weight,  m  the  slope,  and  b  the  intercept.  Thus, 
as  M^ooand  Tg  -»-Tgoo,  the  Tg  at  infinite  molecular  weight,  or  complete  cure,  Tgoofnustbe 
b.  Because: 


where  d  is  the  density  and  c  the  concentration. 


^  "'"g  "'"goo 


^  ■  m  ^''’9  “  ’’’go 


Initially: 


*^0  -  Tg^ 


“  m  “  ''’900  > 


The  fraction  reacted  is  thus  given  by: 


(Cq-c)  ^  "^go  '  “^g 
^0  Tg^  -  Tg^ 


regardless  of  the  kinetic  order  of  the  reaction. 
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For  first-order  reactions: 


-  =  e-M' 

Cn 


where  ki  is  the  first-order  rate  constant,  and  t  is  the  time  passed  since  the  initiation  of  the 
reaction.  Substituting  Equations  B-1  and  B-2  (or  B-3)  into  Equation  B4: 


"^900 
“  "’"soo 


"Tgo  "  "^goo 

PlQ^3  Qf  In  - -  vs  t  give  straight  lines  passing  through  the  origin  with  ki  as  the  slope. 

L'^g'^gooJ 


For  all  other  reaction  orders: 


1  [,n-l  '  ‘ 


where  kn  is  the  nth  order  rate  constant.  Substituting  Equations  B-1  and  B-2  into  Equation 
B-5  gives: 


knt  = 


d^-l  (n-1)  (Tg^-  (Tg^-  Tg^)"''' 


Plots  of  (Tg  -  Tg)^  vs  t  give  straight  lines  passing  through  the  origin  with  slopes  of 


knd^"''  (n-1) 


and  intercepts  of  (Tg^-  Tg^) 


The  rate  constants,  kp,  are  the  same  as  those  in  Appendix  A.  Thus,  by  deterrnining  the 
(density,  m  can  be  determined  from  a  combination  of  a  rate  constant  obtained  from  heat 

evolution  data  and  the  term  knd^’^  (n-D/m^-l  obtained  from  Tg  data  for  the  same 
heat-treatment  temperature.  Thereafter,  the  molecular  weight  of  a  sample  can  be  found 
from  a  Tg  measurement. 


25 


APPENDIX  C 

USE  OF  KINETIC  PARAMETERS  IN  DETERMINING  EXTENT  OF  CURE 

Disappearance  of  a  reactant  at  concentration  c  is  given  by: 


—  =  c\n°e  Ft<T  +  pt)  ^  (C.-i) 

dt  " 


where  t  is  the  time,  n  the  order  of  reaction,  kp®  the  Arrhenius  or  frequency  factor, (c)  E+ 
the  activation  energy,  R  the  gas  constant,  T  the  initial  temperature,  and  p  the  rate  of 
temperature  change.  For  isothermal,  second  order  reactions.  Equation  C-1  maybe  rearranged 
and  integrated  to  give: 


Equation  C-2  may  be  integrated  stepwise  to  give: 


(c)  It  was  assumed  for  ease  of  computation  that  kp®  was  independent  of  temperature. 

(d)  The  series  in  the  third  term  converges  slowly  and  requires  carrying  a  large  number  of 
significant  digits. 
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As  an  example,  samples  of  as-received,  B-staged  resin  were  heat  treated  in  the  following 
manner:  Heated  from  323  to  423  K  at  p  =  1.25  deg  min  1,  held  at  423  K  for  30  minutes, 
and  cooled  from  423  to  323  K  at  p  =  -1.25  deg  min  1.  It  was  assumed  that  the  reaction 
started  when  the  temperature  reached  Tqq  (433  K). 

A  Simpson's-rule  integration  of  Equation  C-2,  with  At  =  240  seconds  for  the  heating  and 
holding  parts  of  the  cycle,  gives: 


—  -  —  =  ^x  240  [K333  K423  +  4(K338  +  K348  +  ■  -  •)  + 

c  Cq  3 


2(K343  +  K353  +  •••)]+  l<2°  1800  e‘ 


14,800/1.987x423  , 


(C-4) 


where  the  K  js  are  the  values  of  k2 


Og-E+/RT 


at  T. 


Equation  A-8  can  be  rewritten  as: 


(C-5) 


The  extent  of  reaction  after  the  1.25  deg  min“l  heat  and  423  K  hold  was  computed  from 
Equations  C-4  and  C-5  to  be  64.0%.  From  Figure  8,  the  Tg  at  that  point  was  399  K.  It  was 
therefore  decided  that  the  reaction  would  probably  cease  when  403  K  was  reached  on  the 
cooling  portion  of  the  heat-treatment  cycle.  A  third  term  was  added  to  Equation  C-4  to 
account  for  the  reaction  occurring  upon  cooling  to  403  K  giving: 


L  -  J_  =  [K333  +  K423  4(K338  +  K348  +  •  •  •)■'■  2  (K343  +  K353  -^  . . .)] 

c  Co  3 


-t-1800  k2°e 


-14,800/1.987x423 


+  ——  [K423+  K403  +  4(K418+  K408)  +  2  K413] .  (C-6) 

o 

The  percent  of  reaction,  calculated  by  substituting  the  result  of  Equation  C-6  into  Equation 
C-5,  was  68.6%.  A  two-minute  anneal  at  423  K  was  needed  to  eliminate  the  very  sharp  endo¬ 
thermic  spike  found  after  the  heat-treatment  cycle,  and  a  Tg  of  407.0  K  was  then  obtained 
in  the  usual  manner.  The  degree  of  cure  can  be  determined  from  Tg  in  two  ways:  (1 )  from 
Equation  B-3,  using  Tgo  =  333.1  K  and  Tgoo  =  434.4  K,  a  Tg  of  407.0  K  corresponds  to  a 
73.0%  cure;  (2)  from  Figure  8,  a  Tg  of  407.0  K  corresponds  to  a  73.6%  cure.  Both  of  these 
ways  compare  very  well  with  the  68.6%  computed  from  rate-of-heat-evolution  data. 


